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Abstract 

The zero forcing number Z[G) of a graph G is the minimum cardinality of a set S of 
black vertices (whereas vertices in V(G)\S are colored white) such that V(G) is turned 
black after finitely many applications of "the (classical) color change rule": a white 
vertex is converted to a black vertex if it is the only white neighbor of a black vertex. 
Zero forcing number was introduced and used to bound the minimum rank of graphs 
by the "AIM Minimum Rank - Special Graphs Work Group". We introduce here a 
probabilistic color change rule (peer) which is a natural generalization of the classical 
color change rule. We introduce a theory of probabilistic zero forcing arising out of the 
peer; the theory yields a quantity Pa{G), which can be viewed as the probability that 
a graph G with an initial black set A will be converted entirely to the color black. We 
also interpret the evolution of the sample spaces of this theory as a Markov process. 
We end with a few basic examples illustrating this theory. 

Key Words: zero forcing, probabilistic zero forcing, probabilistic color change rule, Markov process, 
absorbing chains 
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1 Introduction 

Let G = (V(G), E(G)) be a finite, simple, connected, and undirected graph of order |V(G)| = n > 2. 
For a vertex v € V(G), the open neighborhood of v is the set Ng(v) = {u \ uv € E(G)}, the closed 
neighborhood of v is the set A^H = Ng(v) U {«}, and the degree of a vertex v is deg G (u) = |ATg(u)|; 
we drop G if it is clear in context. 



Let each vertex of a graph G be given either the color black or the color white. Denote by Z the 
initial set of black vertices of G. The "color change rule" changes the color of a vertex v from white 
to black if v is the only white neighbor of a black vertex u; in this case, we say that u forces v and 
write u — > v. The initial black set Z is said to be "zero forcing" if, after finitely many applications of 
the color change rule, all vertices of G are forced to black. The "zero forcing number of G", Z(G), 
is defined as the minimum of the cardinalities of all zero forcing sets of G. 

The notion of a zero forcing set, as well as the associated zero forcing number, of a graph was 
introduced by the "AIM Minimum Rank - Special Graphs Group" (see [I]) to bound the minimum 
rank for numerous families of graphs. Zero forcing parameters were then studied and applied to 
the minimum rank problem in [21 El LT1 El EU H2J H3]- More recently, viewing zero forcing as a 
dynamical process unfolding on a graph, the four authors of [5] and, independently, Hogben et al 
in [TT] have studied the number of steps (or units of time) it takes for an initial zero forcing set to 
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force all vertices of a graph to black. It is also noteworthy that, independent of the forgoing efforts 
and in connection with the control of quantum systems, this parameter has been studied by some 
physicists who called it "the graph infection number" [4J [5j [18] . On the other hand, probability has 
been playing a central role in graph theory ever since the seminal paper by Erdos and Rcnyi [8] . In 
the 2005 Nature article entitled "Evolutionary dynamics on graphs" |15) . Lieberman, Hauert, and 
Nowak studied the question: "what is the probability that a newly introduced mutant generates a 
lineage that takes over the whole population?" 

Let us henceforth use the term "classical zero forcing" to refer to the notion of zero forcing already 
defined. Taking the viewpoint that classical zero forcing is a kind of "graph infection" or "spread 
of a mutation" , it's then natural - indeed compelling - to formulate a notion of "probabilistic zero 
forcing" wherein classical zero forcing is obtained as a special case. To this end, let Z be an initial 
black set of a graph G; here Z ^ 0, but Z is not necessarily zero forcing in the classical sense. With 
respect to Z, we define a "probabilistic color change rule" (hereafter also referred to as "peer") 
in which a black vertex forces each of its white neighbors independently to the color black with 
probability given as follows. Let F(u — ¥ v) denote the probability of the event u — > v (which we 
think of as the "forcing probability of u on u") that the vertex u forces the color of a white vertex 
v to black. We adopt the following 



Definition 1.1. 



]fu<£ZoTv£ N{u) or N[u] C Z 

F(u^v) = l i^nzi otherwise . 

I deg(u) 



We will assume that a vertex u forces each of its white neighbors independently; i.e., we think of 
the action of a black vertex on each of its m white neighbors as a binomial (or Bernoulli) experi- 
ment with probability of success (where by a success we mean the conversion of one of the m white 
neighbors to black) given by the formula above. Since F(u — > v\) = F(u V2) when v\ and 1)2 are 
both white neighbors of u, we will simply write F(u) to specify the forcing probability of u when 
ambiguity is not a concern. Let P{— > v) denote the probability ("conversion probability") that a 
vertex v will be black after one global application (iteration) of peer, given a set of black vertices 
Z C V{G). We thus define P(-» v) = 1 if v G Z and P(-» v) = F{U ueN(v) u -> v) if v G V{G)\Z, 
where F(U ueN ( v) u -> v) = F(m -> v)+F(U ueN ( v )\{ Ul }U —} v)-F(u 1 -> v) ■ F(U ueN ( v )\{ Ul }U -» v); 
i.e., we assume that each white vertex is being acted upon independently by its black neighbors. 

Observe that classical zero forcing where a black vertex u forces its sole white neighbor v is recovered 
as a pair, a black vertex u and a white vertex v, with forcing probability F(u — > v) = 1 and conversion 
probability P(— > v) = 1. From the perspective of modeling a "real world" situation, our formulation 
is sensible: it only takes one affected entity to spread a condition to an entire population. We want 
to mention here that Curtis Miller has independently formulated a number of probabilistic variants 
of zero forcing [16j . But our present definition, as well as our approach and emphasis, is distinct from 
the ones he considered (personal communication). Notice that, once a classical zero forcing (black) 
set Z\ containing an initial black set Zq of G is attained through probabilistic zero forcing, classical 
zero forcing will be in force as a part of probabilistic zero forcing; before then, some combination of 
classical and probabilistic zero forcing occurs. In the rest of this paper, we will at times abbreviate 
probabilistic zero forcing as "PZF" . 



2 Interpretations of PZF 

Let's call a graph G, with vertices labeled 1 through n, a colored graph if each of its vertices receives 
the color black or the color white, and let's regard as distinct two colored graphs G\ and G2 with 
the same underlying graph G when there is a vertex v £ {1, . . . , n} which has one color in G\ but 
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has the other color in G 2 . For a colored graph G, denote by D(G) the set of all colored graphs which 
can be obtained with positive probability from G in one global application of peer. Let a graph G* , 
with initial black set Z, be given. Then PZF generates a countable collection of probability sample 
spaces <£(Go) = {S : k G N}, one sample space S at each step k which can be inductively defined 
as follows: S° = {G }, where G is G* with Z C V(G*) colored black and V(G*) \ Z colored white; 
once S k ~ 1 is defined, S k = (JceS*- 1 D{G'). Let be the probability function defined on S k 
which is derived from the peer. It follows that J2ces k P^ k '(G) = 1 f° r each k G N. The €(Go) just 
defined represents a discrete dynamical system on the underlying graph G, and its evolution is the 
central problem of the PZF theory. 

Let W C V(G*) and Sl^ = {G £ S k : each v G W is a black vertex of G}; if W = {v}, we'll write 
S k . Define P^(W) = f^ces" P {k) (G)- If W u W 2 C V(G*), then it follows that P( k \(-^ Wi) U (-> 

W 2 )) = P (fe) (-> W^i)+P (fc) (-> W 2 )-P^(-> (W1UW2)), since —¥ (W1UW2) = (-)• Wi) n (-> W 2 ). 
It's clear that, relative to an initial black set Z, pW(— > v) = P(—> v), where P(— > w) is previously 
defined. Also notice that p( fc+1 )(— >■ v) > p( fc )(— > w ) by co/or change asymmetry; namely, if a vertex 
v is black in a colored graph G, then w is black for each member of D(G). Suppose for some ko, 
there is no G £ S k ° such that G contains a black vertex v which forces with probability one (i.e., v 
has a lone white neighbor), then S k C S k+1 (as sets) for fc > feo- Since there are at most 2^ V ^- G " 
colored graphs on the underlying graph G*, the sequence S k is eventually constant; i.e., there is 
some M G IN such that S k+1 = S k for any k > M. Note that this docs not mean that p( fc+1 ) = p( fc ) 
for k > M. 

The evolution of our discrete dynamical system is a Markov process. Here is the gross view: Let 
Gi, G 2 , . . . , G 2 " be an ordered listing of the 2" possible colored graphs (possible states of the system) 
undcrlied by a graph G of order n: each G; can be identified with the positive integer i. For each 
k G N, we associate a probability vector of length 2™ whose z-th entry denotes the probability 
that our system is in the state i at time k. The stationary transition matrix A for the system has 
size 2™ x2™, and each of its entries Aij provides the probability that a system in state j at time k 
will transition into state i at time k + 1; Ai t j is thus determined by the peer applied to the state j. 

We observe that, starting with a colored graph Go of order n in S°, by the fact that the black set is 
expanding (or non-contracting) from any G G S k ~ 1 to each G' G D(G), the number of colored graphs 
(and hence the size of the transition matrix) requiring consideration is typically and tremendously 
reduced from 2™. Recall that, in the theory of Markov processes, an absorbing state is one that 
can not be exited once entered; i.e., it is a state i such that A^j = 1; a state i is transient (also 
called non-recurrent) if there is a state j reachable from i, but state i is not reachable from j. If 
there is only one absorbing state jo, then lim A\ • = 1 for each state i (see p. 247 and p. 241 of [I7| . 

fe— >oo " 

also see pp. 58-66 of [14j for proofs). Since our dynamical system has transient states except for one 
absorbing state (the colored graph G& where all vertices of Gb are black), we have the following 

Proposition 2.1. For any colored graph G , lim P' fc '(G&) = 1, where Gb G S k G £(G ). 

k— ¥00 

That is, a colored graph will turn entirely black with probability one as the number of steps k goes 
to infinity. The following is an instructive example. 

Example 1. Consider the complete bipartite graph Ki iTn , where m > 1, which we'll also call the 
star with m pendants; let the center vq be the only black vertex of Go- We have 5° = {Go} and 
S k = {Go, . . . , G 2 m_i} for k > 1; i.e., all possible colored graphs (or states) occur in every sample 
space S for k > 1. Let G& denote the state where all vertices of the star are black. Then Proposi- 
tion l2.1l savs that lim^oo P^ k \Gb) = 1, which can also be seen via a direct, inductive argument on 
the number of pendants m via the following observation: the center vq is the only vertex capable of 
forcing; each G' G D(G) C S k+1 is obtained from G G S k through a binomial experiment wherein 
each white pendant is turned black by the forcing strength of vq; the average forcing strength of vq 
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is monotonically increasing toward 1 as k goes to infinity. 



□ 



In order that probabilistic zero forcing be a natural extension of classical zero forcing, there ought 
to be a probability Pa{G) assigned to each pair (G, A) which indicates the likelihood that the graph 
G will be forced entirely black given the initial black set A C V(G) under finitely many applications 
of the peer; moreover, this probability should satisfy the following conditions: 1) P${G) = 0, 2) 
P A (G) = 1 if and only A is a zero forcing set of G, 3) P A {G) < P B (G) for A C B C V(G). Given a 
pair (G, A) and the associated C(Go), define T k to be the subset of S k such that G 1 G T k if and only 
if the black set of G' contains a classical zero forcing set of G; note that T k ^ for some k G IN by 
the definition of peer. It's clear that the following definition of Pa(G) meets the three requirements. 

Definition 2.2. Let P A {G) = P k «(T k °), where k is the least k such that T k ^ 0. We think of 
P A (G) as the probability that a graph G with an initial black set A C V(G) "will turn entirely black 
in finite time" . 

Example 2. Let G = C n , n > 3, with initial black set A = {vi} C F(G) be given; also let E(G) = 
{viV2, V2V3 . . . ,v n -iv n ,v n vi}. Recall two adjacent vertices form a classical zero forcing set for G. 
Then, in going from S° to S 1 , the only forcing action is given by F(v\ — > v 2 ) — \ = F[v\ —tv n ). 
Only one (call it G w ) of the four colored graphs of S 1 does not contain a classical zero forcing set, 
and thus P A (G) = 1 — P(G W ) = 1 — ^ • | = |. Note the equivalence, trivially here, of P A (G) to 
«a) U (-> «„)) = P«H «2) + «n) " ^ (1) (-> {V2,v n }). □ 




Figure 1: G with initial black set {uq} 



Example 3. Let G be the tree given in Fig[I] with initial black set A = {v n }. The forcing action 
propagates from vq to the three pendants (vertices of degree 1). In Fig[I] the number next to each 
edge uv is F(u —> v)] notice that forcing can only occur in one direction for each edge in this pair 
(G, A). For example, the number | is next to the edge U11U21 because F{v\\ — > V21) = | given that 
the vertex v\\ is already (forced to) black. Likewise, the number next to each vertex indicates the 
probability that each vertex will be forced black given the initial black vertex vq. 

We determine P A (G) by starting with the observation that a colored graph containing a classical zero 
forcing set first occurs in S 2 , since a black set B containing vq is classically zero forcing if and only if 
B n {t>2i, V22} 7^ 0- The top two rows of Figure [2] show S° and S 1 , whereas the third row only shows 
T 2 C S 2 , namely the colored graphs of S 2 which contain classical zero forcing sets. From left to 
right, the first three graphs of the third row form D(G2) and last three graphs of the third row form 
D{Gi), where G2 and G4 are the second and the fourth graphs, respectively, of the second row. Now, 
it's clear that P«(G 2 ) = \ = P (1 )(G 4 ). Thus, P A (G) = P< 2 >(U G 'eT' G ') = ^Kll + ll + ll)] = I 

We can also deduce P A (G) in this case without explicitly enumerating elements of the sample spaces 
£(Go) — {S k : k G N} in the following way. First, still notice that a black set B containing vo is clas- 
sically zero forcing if and only if B n {v 2 i, V22} ^ 0; also notice that P^'(— > v 2 i) = = P^(— > V22), 
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but P< 2 )(^ v 2 i) ^ and p( 2 '(^ v 22 ) ^ 0. Thus, the probability that G will be turned entirely 
black by {vq} under PZF is given by P( 2 )((— > v 2 i) U (— > v 22 j), which must be the same as the 
probability P(E 2 ) of the event E 2 that Dn — > v 2 i or vn — S- v 22 . Let E\ denote the event that 
vq — > vu. Observe Ei C E 2 ; i.e., if E 2 happens, then Ei already happened. Thus, we have 
P{E 2 ) = P{E 2 C\Ei) = P(E 2 \Ei)P(Ei) = (F(vn^v 2 i) + F(vn^v 22 ) - F(v 11 -^v 2 i)-F(v 11 -> 

v 21 )).F(v ^ Vll ) = d + i- 1 •§)•§ = !■ □ 








£\ £\ £\ 

Figure 2: S°, S 1 , and T 2 C S* 2 associated with the pair (G,A) 






The following is a basic question in our PZF theory at this juncture - one which we'll solve for the 
graph depicted in Figure [T] in the case the initial black set is a singleton. 



Question 1. Let P^)(G) = max({PA(G) : \A\ = j}), where A is an initial black set of G and j is 
an integer between 1 and Z(G). What is P(j)(G) for each j £ {1, . . . , Z(G)} and what are the initial 
black sets A attaining P^ (G) for each j? 



Example 4. Let's answer Question [T] for j = 1 when G is the star A"i m with m pendants, and let 
m > 3 for simplicity. Label the center of the star vq and the pendants v± through v m . Given the 
symmetry group of the star, it's enough to consider two cases: 1) A = {uo}, 2) A = {vi}. In case 

1) , F(v ) = i and T 1 ± 0; thus, P {M (G) = P^T 1 ) = ( m ™J ( J-)™"^) = In case 

2) , we have F(v\ — > vq) = 1 and T 1 = since m > 3 by hypothesis. In passing from S 1 to S 2 , we 

have P(« ) = I and T 2 ^ 0; thus, P {vi} (G) = P( 2 )(P 2 ) = (^ZD&'H^) = 
Obviously, P^(G) = P{, Ul j(G) > P{ Vo y(G) when m > 3. We note that it's a straightforward matter 
to determine Pq) (G) for each 1 < j < Z(G) = m — 1 since, for each A of cardinality j, we only need 
to consider the case vq & A and the case Vq £ A. □ 



Example 5. Returning to the graph depicted in FigurcQ] let's consider all cases A4 where |j4,| = 1. 
There are four distinct cases to consider: 1) Ai = {vo}, 2) A 2 = {^23}, 3) A3 = {i>ii}, 4) A4 = {t>2i}- 
Having found in case 1) that Pa 1 (G) = | in the previous example, let's consider case 2). By the 
peer, we have P(— > v 23 ) = P^\— > vi 2 ) = P (2 )(— > vq) = p( 3 )(-> u u ) = 1, and we have in S 4 
that F{v n -> V21) = I = F(v u -> 1-22); we get P A2 (G) = p( 4 )(T 4 ) = P< 4 )((^ v 2 i) U (-»• «»)) = 
§-h§ — §§ = §- In case 3 ): Tl + and ^4 3 (G) = P^T 1 ) = 3 • (| • § • |) + \ ■ \ ■ \ = Zf. In case 
4), since F(v 2 i — > i>ii) =1, S 1 contains only the colored graph where the only black vertices are i>2i 
and vn, and T 1 = 0. Since T 2 ^ 0, P Ai {G) = P< 2 )(T 2 ) = 1 - | • | = |. Therefore, we conclude 
that P ( i)(G) = max({P A (G) : \A\ = 1}) = |. □ 
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